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Abstract 

In this article, we identify the necessary and sufficient conditions for 
the existence of a random field solution for some linear s.p.d.e.'s of 
parabolic and hyperbolic type. These equations rely on a spatial op- 
erator C given by the L 2 -generator of a ci-dimensional Levy process 
X = (Xt)t>o, and are driven by a spatially-homogeneous Gaussian 
noise, which is fractional in time with Hurst index H > 1/2. As an 
application, we consider the case when X is a /3-stable process, with 
(3 € (0, 2]. In the parabolic case, we develop a connection with the po- 
tential theory of the Markov process X (defined as the symmetrization 
of X), and we show that the existence of the solution is related to the 
existence of a "weighted" intersection local time of two independent 
copies of X. 
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1 Introduction 



In 1944, in his seminal article [19], Ito introduced the stochastic integral with 
respect to the Brownian motion, which turned out to be one of the most 
fruitful ideas in mathematics in the 20th century. This lead to the theory of 
diffusions (whose origins can be traced back to [20]), and the development of 
the stochastic calculus with respect to martingales (initiated in [25J). These 
ideas have grown into a solid branch of probability theory called stochastic 
analysis, which includes the study of stochastic partial differential equations 
(s.p.d.e.'s) 

Traditionally, there have been several approaches for the study of s.p.d.e.'s. 
The most important are: the Walsh approach which relies on stochastic in- 
tegrals with respect to martingale-measures (see [33]), the Da Prato and 
Zabczyk approach which uses stochastic integrals with respect to Hilbert- 
space- valued Wiener processes (see [10]), and the Krylov approach which uses 
the concept of function-space- valued solution (see [21]). These approaches 
have been developed at the same time, and nowadays a lot of effort is dedi- 
cated to unify them (see the recent survey [8] and the references therein). 

The fractional Brownian motion (fBm) was introduced by Kolmogorov 
in [23], who called it the 'Wiener spiral", and is defined as a zero-mean 
Gaussian process (B t ) t > with covariance: 

R H (t, s) = E{B t B s ) = l -(t 2H + s 2H - |t - s\ 2H ). 

The parameter H lies in (0, 1), and is called the Hurst index (due to [IE]). The 
case H = 1/2 corresponds to the Brownian motion, whereas the cases H > 
1/2 and if < 1/2 have many contrasting properties and cannot be handled 
simultaneously. The representation of the fBm as a stochastic integral with 
respect to the Brownian motion on M was obtained as early as 1968 (see |26j). 
but the fBm began to be used intensively in stochastic analysis only in the 
late 1990 's. It is the flexibility which stems from the choice of the parameter 
H that makes the fBm a much more attractive model for the noise than the 
Brownian motion (and its infinite-dimensional counterparts). 

Among the fBm's remarkable properties is the fact that it is not a semi- 
martingale. Consequently, Ito calculus cannot be used in this case. A 
stochastic calculus with respect to the fBm was developed for the first time in 
[TT] . Subsequent important contributions were made in [I], [5] and [12]. The 
stochastic integral used by these authors is an extension of the Ito integral 
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introduced by Hitsuda in JT5] (and refined in [2T] and [SI]), and coincides 
with the divergence operator. These techniques are based on Malliavin cal- 
culus. Alternative methods for defining a stochastic integral with respect to 
the fBm exploit the Holder continuity property of its sample paths and are 
based on generalized Stieltjes integrals. We refer the reader to Chapter 6 of 
the monograph jJH] for more details. 

In the present article, we consider the parabolic Cauchy problem 

— (t,x) = Cu(t,x) + W{t,x), t>0,xeR d (1) 
u(0,x) = 0, xeR d , 

and the hyperbolic Cauchy problem 

Cu(t,x) + W(t,x), t>0,xeR d (2) 

o, x eR d 
0, x e R d , 

where £ is a "spatial operator" (i.e. it acts only on the x variable) given by 
the L 2 (IR d )-generator of a (^-dimensional Levy process X = (X t )t>o, and W 
is a Gaussian noise whose covariance is written formally as: 

E[W{t, x)W(s, y)} = \t- s\ 2H ~ 2 f(x - y), 

for some index H > 1/2 and some kernel / (to be defined below) . 

The rigorous definition of the noise W is given in Section [2J At this point, 
we should just mention that the covariance structure of the noise has two 
components: a spatially-homogeneous component specified by the kernel / 
(the example that we have in mind being the Riesz kernel f(x) = Ixl - ^ - "-*, 
with < a < d), and a temporal component inherited from the fBm. This 
becomes clear once we realize that if H > 1/2, Rn{t, s) can be written as: 

R,H(t,s) = aH / / | it — v \ 2H ~ 2 dudv , with olh = H(2H—1). 
Jo Jo 

(The case H < 1/2 has to be treated differently and is not discussed here.) 

The solution to problem (or (J2J)) is understood in the mild-sense, and 
one of the goals of the present article is to give a necessary and sufficient 
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condition for the existence of this solution, in terms of the parameters (H, f) 
of the noise, and the spatial operator C. A similar problem has been consid- 
ered in [H] and [15] in the case H = 1/2 (which corresponds to the white 
noise in time). This motivated us to examine the case H > 1/2. 

The case of the hyperbolic equation with spatial operator C = — (— A) - ' 3 / 2 , 
> 0, driven by a white noise in time was examined in [7] and [9J. In fact, 
these authors consider the much more difficult case of the non-linear equation 
d tt u = Cu + a(u)W + b(u) with arbitrary initial conditions, and Lipschitz 
continuous functions a and b. For the linear equation, it turns out that the 
necessary and sufficient condition for the existence of the solution is: 

where the measure \x is the inverse Fourier transform of / in 5'(IR d ). 

In the case of equations (CQ) and (EJ) driven by a space-time white noise 
(i.e. H = 1/2 and / = So), the authors of [13] have shown that the necessary 
and sufficient condition for the existence of a random field solution is: 

/ — r^n(d£) < oo, (4) 

where is the characteristic exponent of the underlying Levy process A. 
An important observation of [H] is that condition (J4]) can be extrapolated 
in a different context, being the condition which guarantees the existence 
of a local time J* So(X s )ds of the symmetrization X of X. This line of 
investigation was continued in [15] in the case of the parabolic equation (CD) 
with white noise in time, but covariance kernel / in space. Surprisingly, it 
is shown there that condition (j3J) is related not only to the existence of the 
"occupation" local time L t (f) = J * f(X s )ds, but also to the potential theory 
of the process A, when viewed Markov process. 

In the present article, we carry out a similar program in the case of the 
fractional noise in time. More precisely, after the introduction of some back- 
ground material in Section [21 the article is split between the two problems: 
Section [3] is dedicated to the parabolic problem (pQ), while Section 0] treats 
the hyperbolic problem (j2J). For the parabolic problem, we discuss three 
things: (i) the existence of a random field solution (Section 3.1); (ii) the 
maximal principle which gives the connection with the potential theory of 
Markov processes (Section 3.2); (iii) the relationship with the "weighted" 
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intersection local time L tt n(f), defined by 



Lt,H(f) = a H 



L 



r — 



s\ 2H - 2 f{Xl-Xl)drds, 



where X 1 and X 2 are two independent copies of X (Section 3.3). For the 
hyperbolic problem, we only discuss the existence of a random field solution 
in the case when is real- valued (i.e. X is symmetric). 

Unlike the case of the white noise in time, it turns out that for the frac- 
tional noise, the conditions for the existence of the solution are different for 
the parabolic and hyperbolic problems. These conditions are: 



in the hyperbolic case. This phenomenon was observed for the first time in 
[I] for the wave and heat equations. As an application, we discuss the case 
when X is a /3-stable process with (3 G (0,2], and hence = cp\£\ p . In 
this case, conditions flH]) and ([6]) turn out to be generalizations of ([3]). 

The fact that the fractional noise induces a connection with the weighted 
intersection local time was also used in [T7], in the case when / = 5q. In 
[3], it is shown that the existence of the exponential moment of L tt jj(f) is 
closely related to the existence of the (mild) solution of the heat equation 
with multiplicative noise. 

We now introduce the notation used in the present article. The Fourier 
transform of a function if G L 1 (M d ) is defined by: 



It is known that the Fourier transform can be extended to L 2 (M. d ) (see e.g. 
[13]). Plancherel theorem says that for any (p, ip G L 2 (M. d ), 





in the parabolic case, respectively, 
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Let S(R d ) be the Schwartz space of rapidly decreasing infinitely differ- 
ential) le functions on M. d . A continuous linear functional on is called 
a tempered distribution. Let S'(M. d ) be the space of tempered distributions. 
The Fourier transform FS of a functional S G S'(M. d ) is defined by: 

(FS, <p) = (S, T(p) , Vy> G S(R d ) . 

We refer the reader to [30] for more details about the space iS'(IR d ). 

2 Preliminaries 

In this section, we introduce some background material about the Gaussian 
noise W and the Levy process X. 

2.1 The Gaussian noise 

As in [6], we let / : M. d — > [0, oo] be a measurable locally integrable function 
(or a kernel). We assume that / is the Fourier transform in £'(1^) of a 
tempered measure /i, i.e. 

/ f(xMx)dx= [ JV(0MO, V^5(l d ). (7) 

It follows that for any (p, ip G 5(lR d ), 

/ / y{x)iP{y)f{x-y)dxdy= [ (8) 

JR d JR d JR d 

Similarly to [I], we let £ be the set of elementary functions of the form 

h(t, x) = (p(t)^{x), t>0, x G M. d , 

where is a linear combination of indicator functions l[o, a ] with a > 0, and 
i) G S(R d ). We endow £ with the inner product: 



(hi, h 2 )nr — 




2 f(x — y)h\(t, x)h2(s, y)dxdydrds. 



Let 7-LV be the completion of £ with respect to the inner product (•, -)uv- We 
note that the space T-LV may contain distributions in both t and x variables. 
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We consider a zero-mean Gaussian process {W(h); h G £} with covariance 
E(W(h)W(h 2 )) = (hxMuv- 

The map H- ^(/j) is an isometry between £ and the Gaussian space 
of W, which can be extended to HP. This extension defines an isonormal 
Gaussian process W = {W(h); h G HP}. We write 

poo i* 

W(h)= / h{t,x)W{dt,dx), for any he HP. 

JO JR d 

This defines the stochastic integral of an element h G HP with respect to 
the noise W. 

2.2 The Levy process 

As in [15], we let X = (X t )t>o be a <i-dimensional Levy process with charac- 
teristic exponent Hence, Xq = and for any t > 0, 

E( e -*' Xt ) = e- m ^\ for all £ G M d . (9) 

By the Levy-Khintchine formula (see e.g Theorem 8.1 of [29|), 

^(0=^7-^ + ^- I (e-^ x -l + it-xl {lxl<1} )v(dx), 

where (7, A, v) is the generating triplet of X. We observe that for any £ G M. d , 

Re* (f) = £ T A£ + / [1 - cos(£ • > 0, and 

Im#(£) = 7-f+ / [sin(£ • x) -£ ■ a;l { | x |<i } ]z/((ix). 

X is a homogenous Markov process with transition probabilities: 

Q t (x; 5) = P(X S+ , G P|X S = a;) = P(X t eB-x), 

for any x G M d and Borel set B C R d . Let (P f )oo be the associated semi- 
group, defined by: 

(P t <f>)(x) = [ <f>(y)Q t (x; dy) = E[<f>(x + X t )}, 



for any bounded (or non-negative) measurable function : M. d — > M. 
Let C be the L 2 (M d )-generator of (Pt)t>o, denned by: 

£0 = li m ^tlA i n J} md\ (if it exists) . 
t-to t 

Note that £0 exists if and only if G Dom(£), where 

Dom(£) = {0 G L 2 (R d ); (J^0)^ G L 2 (R d )} 

(see p. 16 of [15J). Moreover, £ can be viewed as a convolution operator with 
Fourier multiplier TL = i.e. for any G Dom(£) and £ G M. d , 

Hem) = -*(0^(0- 



3 The Parabolic Equation 

In this section, we assume that the law of X t has a density denoted by p t . This 
assumption allows us to identify the fundamental solution of d t u — Cu = 0. 

To see this, note that Qt(x; •) has density p t (- — x), and P t (f> = 4> * Pu 
where pt(x) = p t (—x). Since the solution of the Kolmogorov's equation 
d t u(t,x) = £u(t,x) with initial condition u(0,x) = u (x) is 

u{t,x) = (P t u )(x) = / u (y)p t (y - x)dy, 

JM. d 

it follows that the fundamental solution of dtu — Cu = is the function: 

G{t,x) = Pt (-x), t>0,xeR d . 
From (jUJ), we obtain that: 

J-G(t, -)(0 = / e^ x Pt (x)dx = E(e^ Xt ) = e"'*^. (10) 



3.1 Existence of the Random-Field Solution 

There are two equivalent ways of defining a random field solution for problem 
(CQ). Similarly to [4J, one can say that the process {u(t,x);t > 0,x G R d } 
defined by: 

u(t,x)= / / G(t — s, x — y)W(ds, dy) 

J0 JR d 
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is a random field solution of ([I]), provided that the stochastic integral above 
is well-defined, i.e. the integrand 



.XJ 



3 (s,y) H> g t x(s,y) = l[ ,t](s)G(t-s,x-y) belongs to UV. (11) 



Since gt x satisfies conditions (i)-(iii) of Theorem 2.1 of [3], to check that 
g tx £ UV, it suffices to prove that: 

I t -=a H / / Tg tx (r, ^(O^Qtxis, -)(0l r ~ s\ 2H ' 2 drds/i(d^) < oo. 

JR d JO JO 

In this case, 

E\u(t,x)\ 2 = E\W(g tx )\ 2 = \\g tx \\ 2 nv = h- 
Note that by (HUD, I t = J Rd N t (£)fi(d£), where 

N t (0 = a H f f e- r ^e~ s ^\r - s\ 2H ' 2 drds. 
Jo Jo 

Therefore, the question about the existence of a random field solution of (pQ) 
reduces to finding suitable upper and lower bounds for N t (£). 

Alternatively, the authors of [2] suggest a different method for defining 
the random field solution of ([1]), which has the advantage that can be ap- 
plied also to the hyperbolic problem (T5]) (for which one cannot identify the 
fundamental solution G). Since this is the method that we use in the present 
article, we explain it below. 

We say that {u(t, (f)',t>0,<fE S{R d )} is a the weak solution of ([T]) if 



u(t,(p)= / / I / G(t — s, x — y)(p(x)dx J 

JO JR d \JR d J 



G{t — s, x — y)ip{x)dxj W(ds, dy). 

Note that the stochastic integral above is well-defined (as a random variable 
in L 2 (fl)) if and only if the integrand 

t+xK d 9 (s,y) \-> h tt9 (s,y) = l[a,t](s)(<p *Pt-s)(y) belongs to UV. 

Since h tlip satisfies the conditions (i)-(iii) of Theorem 2.1 of [4J, to check 
that h tjtp G UV it suffices to show that 

p pco poo 

h, v '=oi H / / \r- s\ 2H ~ 2 Tht^{r, ^(QFht^s, -)(£)drdsiJ,(d£) < oo. 

jR d Jo Jo 
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In this case, we have 

E\u(t^)\ 2 = E\W{h t ^)\ 2 = \\h t J\ 2 HV = I t>ip . 
Since both (p and pt-s are in 

T(<P*Pt-s)(0 = ^(O^Pt-(0 =^(0e- (t - a) * ( «, 

and hence 

h v = [ N t (0\^(0\ 2 fi(dO- 

Using the trivial bound iV f (£) < t 2H (which is obtained using the fact 
that | e -'*®| = e _dto *® < 1 for all s > and £ G we get: 

it lV < t 2H [ I^(0IV(^0 < 00 for a11 e S(M d ). 

This shows that u(t, (p) is a well-defined random variable in L 2 (Q). 

We continue to explain the method of [H]. We endow iS(R d ) with the 
inner product: 

(ip,ip)t = E(u(t,ip)u(t,ip)). 
We denote by || • \\t the norm induced by the inner product (*, -)t, i.e. 

M 2 = E\u(t^)\ 2 = [ N t (0\^(0\ 2 KdO- (12) 

JRd, 

Let M t be the completion of iS(M d ) with respect to (•, -) t and M = n t> oM t . 
We say that (JX^ has a random field solution if and only if 

S x G M for all x G M d . (13) 

The random field solution is defined by {u(t, x) = u(t,6 x );t > 0,x G 

To prove ([TBI , we introduce the space 2 = r\ t> o2 t , where Z t is the 
completion of «S(M d ) with respect to the inner product [-,-]t defined by: 

^ * = I ( l/ t + L({) ) " ^«>»< d « = : *>■ « 
We denote by ||| • the norm induced by the inner product [■, i.e. 
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By Lemma [A.2I (Appendix A), for any s,t > 0, there exist some positive 
constants Ci(s,t) and c 2 (s,t) such that for any <p G <S(IR d ), 

Cl (s, t) 2H £{s- <p) < £{t- <p) < c 2 (s, t) 2H £{s- <p). 

Therefore, the norms ||| • \\\ t and ||| • ||| s are equivalent and Z t = Z s = Z. 

The idea for proving (fT3l) is to show that any norm || ■ || t is equivalent to a 
norm ||| • \\\ p (t), for a certain bijective function p : M + — > R + . From this, one 
infers that M t = Z p m for any t > 0, and hence M = Z = Z\. Condition f fT3|) 
becomes 5 X G Z\ for all x G M d , for which one can find a natural necessary 
and sufficient condition (see Theorem 13 .31 below) . In the case of the parabolic 
problem (PQ), it turns out that p{t) = t. (We will see in Section 0] that for the 
hyperbolic problem p(t) = t 2 .) 

The next theorem is the main result of the present section, and gives 
the desired upper and lower bounds for A^(£). Unfortunately, for the lower 
bound, we had to introduce an additional condition of boundedness on the 
ratio between the imaginary part and the real part of the characteristic ex- 
ponent A similar difficulty has been encountered in [22] for obtaining 
a lower bound for the "sojourn operator". Our condition f|T6l) is similar to 
condition (3.3) of [22], and is trivially satisfied when \1/ is real- valued. 

We use the following inequality: there exists a constant bu > 0, such that 

/ r \ 2H 

\2H-2, 




«h / / |</>(r)||^(s)||r 



2 drds<b 2 H (^J \<p(s)\ l/H ds^j (15) 



for any ip G L l l H (Jsl). This inequality was proved in [2Z] and is a consequence 
of the Littlewood-Hardy inequality. 

For complex-valued functions (p, we define: 



t rt 



MI«(o.t) -= a H / <p(r)<p(s)\r - s\ drds = \\Re<p\\ Hm + \\Im<p\\ nm . 
Jo Jo 



Theorem 3.1 For any t > and £ G 

/ 1 \ 211 



l/t + Re*(0 

where Ch = H 2H b 2 H e 2 . If in addition, there exists a constant K > such 
that: 

|Im#(£)| < #Retf(0, V(el d . (16) 
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pd 



then, for any t > and £ e 

m)>c H , K ( 1/t+ l m ) 2 ", 

where Ch,k is a positive constant depending on H and K . 

Proof: For the upper bound, we note that N t (£) can be written as 

iV t (0 = a H f e - rRc *(O e sRcnO\ r _ s | 2 ^- 2 cos [(r - s)lm^^)]drds. 
Jo Jo 

Using the fact that | cosa;| < 1 and e _rIte *^ < e t / A e _r ( 1 / A+Ba *t € W for any 
r G [0, t], we get: 

iV 4 (0 < e 2t ' x a H f /'* e ^( 1 A+Re*«)) e -.(i/A+iu,*«))| r _ s |2H-2 drtfe 
Jo io 

By (USD, it follows that: 

2H / tj \ 2H 



1/A + Re*(£) 



The conclusion follows by taking A = i. 

For the lower bound, suppose first that tRe\I/(£) < a, for some constant 
a = a K e (0, 1) such that Ka < tt/2. By (fT6|) . 

t|lm*(0| < A"tRe*(0 < £Ta < |. 

Using the fact that e~ x > 1 — x for x > 0, we obtain: for any r e [0, t], 
e -rRe*«) > x _ rRe ^(£) > 1 - a. 

Since cos is decreasing on the interval [0, |], for any < s < r < t, 
cos[(r - s)|Imtf(f)|] > cos[t|Imtf (f)|] > cos(Ka) > 0. 

Therefore, 

N t (£) = 2a H f [\- rRe ^e- sRe ^\r - s) 2H ~ 2 cos[{r - s)\lm^{0\]dsdr 
Jo Jo 

ft rr 

> (1 - a) 2 cos(Ka)2a H / / (r - s) 2H ~ 2 dsdr 

Jo Jo 



1 - a) 2 cos(Ka)t 2M > (1 - a) 2 cos(ATa 
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l/t + Re*(0 



where for the last inequality we used the fact that t > pr^^iy • 
Suppose next that tRe\I/(£) > a. Note that 

iV t (OH|e-* (e) llW 

Using Lemma lB.ll (Appendix B) for expressing the "H(0, t)-norm of the 
exponential function in the spectral domain, we obtain: 

f sin 2 [(T + Im*(fl)t] + {e-»°*«> - c03 [(t + Im^KMP , ,_ (2g _ 1) 
We denote 

T = «!«*(£) and 6=^. 
Using the change of variable r' = r/Re\l/(£), we obtain that: 

*«> = ps^p /.rSJ^w + ^ MI<iT - (17) 

where /r(' r ) = sin[(r + 6)T] and <7t(t) = e _T — cos[(r + 6)T]. 
From the proof of Lemma IB. II (Appendix B), we know that: 

:[fHr)+9Ur)] = \^o,Mr)\ 2 , 



1 + (r + b} 



where tp(x) = e~ x ( 1+lb \ 

Let p > K be positive constant whose value will be specified later. Since 
the integrand of (fTTI) is non-negative, the integral can be bounded below by 
the integral over the region |r| < p. In this region, \r\^^ 2H ^ > p~( 2H ~ l > . 
We obtain: 

s tSp ( /(t) - J M> _ e T^hw imT)+9l{T)]dT ) • (18) 

where 

7(T) := f — -4-— [/|(r)+flg.(r)]dr = 2vr T |e^ 1+lfe )| 2 rfx = 7r(l-e- 2r ), 

1 + ( r + or 

(19) 

by Plancherel's theorem. 
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Using (1191) . we obtain the lower bound: 

J(T) > tt(1 - e- 2a ), since T > a. (20) 

To find an upper bound for the second integral on the right-hand side of 
(1T5|) . we use the fact that: 

/|(t) + (#(t)<5, Vr G K. 

It follows that: 

/• /|(r)+^(r) /• 5 lOp lOp 

i M >P 1 + (r + &) 2 " 7 M > P (r + fe) 2 ^ p 2 - 6 2 " p 2 - J^ 2 ' ^ 

since |6| < K (by (TT6|) ) . We choose p = px large enough such that 
Cjc:=Ir (l- e --)--i^_>0. 
Using (jlgj). PD]1 and (Bill, we obtain: 

^(0 > gic g^TT^ > C K c H p-^-V 



The conclusion follows, letting 

C H)if = min {(1 - a) 2 cos(tfa), C K c H p~ [2H ' l) } . 

□ 

The following result is an immediate consequence of Theorem 13.11 

Corollary 3.2 a) For anyt>Q,tpE S(R d ), 

E\u(t,<p)\ 2 < C H S(t; if ), 

where C H = H 2H b 2 H e 2 . Hence, M t D Z t for all t > 0, and M D Z. 
b) If W holds > then foranyt>0,cpe S(R d ), 

E\u(t, V )\ 2 >C H , K £(t;v), 

where Ch,k is a positive constant depending on H and K. Hence, M t = Z t 
for all t>0, and M = Z. 
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Proof: We use Theorem 13.11 and the definitions (1121) and (1141) of the norms 
|| • \\ t , respectively ||| • \\\ t . □ 

The next result gives the necessary and sufficient condition for 5 X G Zi 
for all xeR d . 

Theorem 3.3 In order that 5 X G Z\ for all x G M. d , it is necessary and 
sufficient that condition |3]) holds. 

Proof: Suppose first that ([5]) holds. To show that 5 X G Z\ for all x G M. d , 
we use an argument similar to the proof of Theorem 2 of [6] . 

Let Z Q be the set Schwartz distributions <p such that J-(p is a function 
and 

iii^iii.-/„XttrW)) 2 ' ,| ^ K)|2 '' ( ^ < ° o - 

Note that S(R d ) C Z and the definition of ||j • || |i agrees on S(R d ) with 
the one given by (JT4l) . Therefore, to show that a distribution (p G Zq is in 
Z\, it suffices to show that there exists a sequence (y? n )n>i C iS(IR d ) such 
that \\\<fin — ¥>|||i 0. We apply this to ip = 5 X . In this case, J-<p(C) = e~^' x , 
= 1 f° r an £ e an d I W^P |||i coincides with the integral of ([5]). 

Let (p n = (p * (fi n E S(M. d ), where (fi n (x) = n d <fi{nx) and (fi G 5(lR d ) is such 
that (fi > and J Rd (fi(x)dx = 1. Then J-'ipniO = J^fiO-^^niO and 

by the Dominated Convergence Theorem, since |J-"0 n (OI — 1 f° r a ^ £ ^ 

For the reverse implication, suppose that 5 X G Z\ for all x G lR d . To show 
that (jSJ holds, one can use the same argument as in the proof of Lemma 4.2 
of [IS] . We omit the details. □ 

The following result concludes our discussion about the existence of a 
random-field solution. 

Theorem 3.4 (Existence of Solution in the Parabolic Case) 

a ) If (GP holds, then equation (QP has a random-field solution. 

b) Suppose that [To]) holds. If (QP has a random-field solution, then |3]) 
holds. 
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Proof: a) Suppose that <^ holds. By Theorem 13. 3} 8 X G Z\ for all x G M. d . 
By Corollary O a), Z x = Z C M. Hence (pSD holds. 

b) Suppose that (JS) holds. By Corollary Ob), M — Z — Z\. Hence 
^ G £1 for all x G M d . By Theorem E3\ © holds. □ 

Example 3.5 (Stable processes) Suppose that C = — (— A)^ 2 for /3 G (0, 2]. 
Then (X t ) t > is a rotation invariant strictly /3-stable process on M. d , and 
#(£) = Cp \(\P ( S ee Theorem 14.14 in [23] and Example 30.6 in [2S]). It can 
be shown that (X t )t>o is subordinate to the Brownian motion on M. d by a 
strictly (/3/2)-stable subordinator (see Example 32.7 of [29]). 
In this case, condition fl5]) becomes: 

I(iw) ! "« <M ' (22) 

We consider two kernels: 

(i) f(x) = Ca^lxl - ^ - ") for < a < d. In this case, [i(d£) = \i\~~ a d£ t (see 
p. 117 of |32j), and condition ( I2"2"j) is equivalent to 



2H(3 >d-a. (23) 

(ii) fix) = nti(«^l^| 2Hl - 2 )- ^ this case, »(d£) = nLi^l&r^-^ 
and condition f )22|) is equivalent to 



d 

2H/3 > d-J2( 2H i- !)■ 

i=i 

Remark 3.6 (Fractional Powers of the Laplacian) As in [7] and [9], we can 
consider also the case C = —(—A)' 3 / 2 for arbitrary j3 > 0, even if there is no 
corresponding Levy process whose generator is C. Note that the fundamental 
solution G of dtu — Cu = exists and satisfies: 

FG(t,0 = exp(-cpt\Zf). 

Using Theorem 2.1 of [3] and estimates similar to those given by Theorem 
13.11 above, one can show that a random field solution of (PQ) (in the sense of 
[3]) exists if and only if (J22D holds. 
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3.2 A Maximum Principle 

Throughout this section, we assume that p t G L 2 {R d ) for all t > 0, and 

/x has a (non-negative) density g>, (24) 

i.e. / is a kernel of positive type (see Definition 5.1 of [22]). 

We consider the symmetric Levy process X = (X t )t>o defined by: 

X t := X t — X t , 

where (X t ) t > is an independent copy of (X t ) t > . We denote by (Pt)t>o the 
semigroup of (X t ) t > , i.e. 

(P t <j>)(x) = cj){y)p t {x -y)dy, 

where pt = Pt * Pt- From ffTO]) . it follows that J r pt($,) = e~ 2 * Re *^. 
Let (-Rq,)q,>o be the resolvent of (Pt)t>o, i-e. 

(R a (f))(x) = / e- QS (P s 0)(x)rf S . 
Jo 

The following maximum principle has been obtained recently in [15] : 

(iW)(0) = sup(R a f)(x) = T(a) := / — -L-^de). (25) 
xgRd a + 2ReW(£) 

Remark 3.7 Recall that / = Tg in <S'(M ). The authors of [15] work with 
the Fourier transform J 7 / instead of g, which introduces an additional factor 
(2ir)~ d . To see this, note that by the Fourier inversion theorem on iS(lR d ), 
relation (J7J) becomes: for any G iS(lR d ), 

/ 0(£M£R = / f{x)jFcj>{x)dx. 

JRd (27r) d J R d 

This shows that # = {2ir)- d J : f in 5'(M d ). 

Note that T(a) < oo for all a > if and only if T(a) < oo for some 
a > 0. An important consequence of ( |25|) (combined with the results of [6]) 
is that the potential-theoretic condition: 

(R a f)(0)<oo for all a>0 
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is necessary and sufficient for the existence of a random field solution of (pQ), 
when the Gaussian noise W is white in time (i.e. H = 1/2). 

In the present article, we develop a maximal principle similar to (j25|) . 
which has a connection with the existence of a random field solution of ([I]), 
when the noise W is fractional in time. 

We define the following "fractional analogue" of the resolvent operator: 

pCsO POO 

(R a , H <P)(x) = a H / \r- s\ 2H - 2 e-^ r+s \P r+s( f>)(x)drds, 
Jo Jo 

and we let 

I* poo poo 

T H (a) :=a H \r - s\ 2H - 2 e^ a+me ^ r+ ^drds^dO. 

</R d Jo Jo 

As in [22], we assume that / satisfies the following condition: 

f(x) < oo if and only if x 7^ 0. (26) 
Under this condition, the following harmonic-analysis result holds: 

/ / ( p(x)^(y)f(x-y)dxdy= f ( 27 ) 

jR d jR d jR d 

for any non-negative functions (p,ip £ ^(W 1 ) (see Lemma 5.6 of [22J). 
Theorem 3.8 (A maximum principle) If j[2o)) holds, then for any a > 0, 
(R a ,Hf)(0) = sup (R a:H f)(x) = Tff(a). 

x£R d 

The proof of Theorem 13.81 follows from Lemma 13.111 and Lemma 13.121 
below. Before this, we need some intermediate results. Let Co(lR d ) be the 
space of continuous functions which vanish at infinity. 

Lemma 3.9 For any <p G 5(IR d ) ; we have: 

a) f*<pe C (R d ) and R a , H (f * <P) G C (R d ) for any a > 0; 

b) f * G L 2 (R d ) and F(f * 0)(£) = (2vr)^0(O^(O- 
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Proof: a) Since / is tempered (i.e. f(x) < C(l + | a: | ) ^ for all x G M. d , for 
some k > 0, C > 0), the function / * is well-defined. By (J7J) and (121]) . 

(f*<f>){x)= f f(y)cf>{x-y)dy= [ e^*JWg)g(Z)d£. (28) 

Since g is tempered, (J-(f))g G L p (K d ) for any p > 1. By Riemann-Lebesgue 
lemma, / * G C (M d ). Finally, we note that R^ H : C (M d ) -> C (R d ). 

b) By (EHD, f*<j> = Fh, where ft, := (7^) 5 G L 2 (R d ). Hence, / * G 
L 2 (M d ). By the Fourier inversion formula in L 2 (M d ) (see e.g. p. 222 of [T3]). 

J-(/ * 0)(£) = (27T) d MO = (27T) d .F0(O<?(O- 

□ 

Lemma 3.10 For any G <S(M d ) and x G M d , 

(Ra, H (f*<f>))( x ) = a H [ e-* x T^) f \r-s\ 2H - 2 e~ ia+2Re ^ mr+s) drdsij(dO. 

Consequently, for any G 5(M d ) ||0||i = 1, 

|(iUff(/*0))(aOI<Tir(a), Vx G M d . 
Proof: Since -P r +<j = -P r P s , we have: 

(P r+s (/*0))(x) = / / {f*<f>){y-z)p r (x-y)p a (z)dydz. (29) 

jR d JR d 

Using Lemma [C II (Appendix C) with (p = f *<p, ipi = p r (x — ■) and ^2 — Psi 

(P r+S (f * 0))(x) = / e-^-^Cl) S) Hf*<i>)l£)<%- 

The result follows using Lemma 13.91 b). the fact that J-p r (£) = e~ 2rRe *^, 
and Fubini's theorem. □ 

Lemma 3.11 For any a > 0, 

Tff(a) = sup(R ajH f)(x) = limsup(R a>H f)(x). 

x€R d 
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Proof: The proof is similar to Proposition 3.5 of [15]. By Fatou's lemma 
and Lemma I3.10[ for any x G M. d , 

(R a ,Hf)(x) < \immf(R a , H (f * <p n ))(x) < T H (a), 

n— >oo 

where (0 n )n>i is a sequence of approximations to the identity, consisting of 
probability density functions in Hence, 

swp(R a , H f)(x) < T H (a). 

xeR d 

For the reverse inequality, we let (j> n (x) = (2ir)~ d / 2 n d / 2 exp(— n\x\ 2 /2). By 
Lemma 13.101 

{R a>H {f*<t>n))(0) = Oi H ! e-l«l 2 /( 2 ™) f \r-s\ 2H ~ 2 e-( a+2Re ^ r+s ^drds^dO, 

JR d JR2_ 

and therefore, by applying the monotone convergence theorem, 

lim(^ aiH (/*0 n ))(O) = T H (a). (30) 

n— too 

Using f )29|) and the symmetry of the function <f) n , we obtain: 
(P r +. a (/*^ n ))(0)= / {P r+S f){x)(j )n {x)dx. 

Therefore, 

(R a Af*^nW) = [ (R a , H f)(x)(l> n (x)dx<sup(R a)H f)(x). (31) 
Jm d xm d 

From fl30|) and fl31|) . we infer that Tjy(a) < sup xgR d(-R Q , i ///)(x). 

The last assertion follows by taking <ft n with the support in the ball of 
radius 1/n and center 0. □ 

Lemma 3.12 /// satisfies $M), then for any a > 0, 

(R a , H f)(0) = T H (a). 
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Proof: Using ( 12"T|) . we have: 

(W)(0) = / / p r {x)p s {y)f{x-y)dxdy= f e~ 2 ^ Ke ^ g^. 

The conclusion follows from the definitions of (^ Q) ^/)(0) and Y#(a). □ 
To investigate the connection with the parabolic problem (pQ), we let 



By Lemma IA.1I (Appendix A), T#(a:) < oo for all a > if and only if 
T^-(a) < oo for some a > 0. 

The following result gives the relationship between Y#(a) and T* H (a). 

Lemma 3.13 For any a > 0, 

c a , H T* H (a) < T H {a) < b 2 H H 2H T* H (a), 
where c a , H = 2^ 2H+2 ^(a A I) /{a + 3/2)} 2H . 



Proof: The second inequality follows by (TT5I) . For the first inequality, we 
note that, since the integrand from the definition of Y#(a) is non-negative, 
the integral drds over [0, oo) 2 can be bounded below by the integral over 
[0, l] 2 . By Proposition 4.3 of [I], for any t > and A > 



^'| r - sr -V^ ds >I( t -A 1 )(0 2ff ( I7 ^) 2H . 



Applying this result for t = 1 and A = a + 2Re\I/(£), we obtain: 

1 r l /I \ 2^+2 / -, \2H 

| r _ s |2H-2 e -(«+2IU ) »«))(r+.) drds > ( L ^ 



o jo 
Hence, 



2 J V 1 / 2 + a + 2Re ^(0 



I N 2H+2 



T H (a)> (-) T^(a + l/2)>c Q ,^(a) 



where we used Lemma IA.1I (Appendix A) for the second inequality. □ 

Recall that by Theorem I3.4[ condition (j3J) is the necessary and sufficient 
for problem to have a random field solution. As a consequence of the 
maximum principle, we obtain the following result. 
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Corollary 3.14 Suppose that f satisfies l[26)) . Then |3]] is equivalent to 

{Ra,H f ){0) <oo for any a > 0. (32) 

Proof: By Theorem 13.81 (I3"2"j) holds if and only if T#(a) < oo for any a > 0. 
By Lemma [3.131 this is equivalent to T* H (a) < oo for any a > 0, which in 
turn, is equivalent to (jSJ) (i.e. T* H {2) < oo), by Lemma lA.ll (Appendix A). 
□ 

3.3 Connection with the Intersection Local Time 

When the noise W is white in time, the authors of [H] and [15] noticed an 
interesting connection between the existence of a random field solution of 
problem ([rj and the existence of the occupation time 

L t (f)= f f(X s )ds. 
Jo 

In this section, we develop a similar connection in the case of the fractional 
noise, by considering the "weighted" intersection local time: 

LtAf) = <*h I f \r - s\ 2H - 2 f(Xl - X*)drds, 
Jo Jo 

where (X^) t >Q and (X^) t > are two independent copies of (X t ) t > . 

Clearly for any fixed t > 0, E[L tt H{f)] < oo is a sufficient condition for 
Lt,H(f) < °° a - s -> but the negligible set depends on t. Our result will show 
that under condition (13"2"|) . L t)H (f) < oo for all t > a.s. To motivate this 
result, we consider first an example, in which we proceed to the calculation 
of E[Lt 7 H(f)) m a particular case. 

Note that X} - X 2 S = X r+S for any r, s G [0, t], and therefore, 

E[L t , H (f)\ = a H f f \r - s\ 2H - 2 E[f(X r+s )}drds. (33) 
Jo Jo 

Example 3.15 (Stable processes) Refer to Example 13.51 Since (X t )t>o is 
self-similar with exponent 1//3 (see Theorem 13.5 of [29]). 

X^Air + sf^X,. 
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Suppose in addition that f(x) = \x\ ( d °^ for < a < d. Then E[f(X r+s )] = 
E\X r+s \-( d - a ^ = c a4 (r + s)-^-")// 3 , where c a>d = ElX^-^. By (155]), 

E[L t>H (f)) = a H c a4 f f \r - s| 2 ^ 2 (r + s^-^drds. 
Jo Jo 

One can see that E[L ty H{f)} < oo for any t > if and only if (|23|) (or 
equivalently, ([5])) holds. But by Corollary 13. 14[ fl5]) is equivalent to ( 132|) . 

The previous example shows that the existence of L t ^ H (f) is related to 
the potential-theoretic condition (13"2"j) . which is in turn the necessary and 
sufficient condition for the existence of a random field solution to problem 
([1]) (by Theorem 13.41 and Corollary I3.14p . We will see below that this is a 
general phenomenon. For this, suppose that 

-^o = X l' x o = X2, 

and let P Xi be the law of X % for i — 1, 2. Then -P^,^ = P Xl x -Px 2 is the law 
of (X^X 2 ). We denote by E XltX2 the expectation under P XuX2 . 

The next result shows the existence of the intersection local time L t) H(f) 
under condition ( 1321) . 

Theorem 3.16 (Connection with the Local Time) Suppose that f satisfies 
j[2h)) . If (Elj] holds, then for any xi,x 2 £ R d , 

P XuX2 (L t:H (f) < oo for allt>0) = l 

( log L tjH (f) \ 
P XUX2 hmsup < =1. 

V t-+oo t J 

Proof: We follow the lines of the proof of Theorem 3.13 of [15J. Since / is 
non-negative, it follows that for any t > 0, 

poo roc 

e- 2at L ttH (f) <a H / e- a{r+s) \r-s\ 2H - 2 f(Xl-X 2 s )drds. (34) 
Jo Jo 

Note that 

E xux 2 [f( x l- X s)]= / / f{y-z)Pr(%i-y)Ps(x2-z)dydz = (P r+s f)(x 1 -x 2 ). 

jR d jR d 
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Taking supremum over t, and expectation with respect to P Xl , X2 m 
obtain: 

E XuX2 [sup(e- 2at L t , H (f))} < (R a ,„f)(x 1 -x 2 ). 
t>o 

From here, using Theorem 13.81 and condition (I3"2l . we infer that: 



sup E XuX2 [sup(e 2at L t)H (f))} < sup(R a , H f)(x) = {R a , H f)(0) < oo. 
The result follows. □ 

4 The Hyperbolic Equation 

In this section we consider the hyperbolic problem (J2J). Throughout this 
section, we assume that X is symmetric, i.e. 

Imtf(0 = 0, for all £ e R d . (35) 

Since Re\l/(£) = ^(O? we use t ne notation to simplify the writing. 

To define the weak solution, we cannot use the same method as in the 
parabolic case, since in general, we may not be able to identify the funda- 
mental solution G of duu — Cu = 0. Note that in some particular cases, we 
are able to identify G (see Remark [4771 below) . 

To circumvent this difficulty, we use the method of [T4] , whose salient 
features we recall briefly below. Consider first the deterministic equation: 

d 2 u 

— (t,x) = jTu(t,x) + F(t,x), (36) 

with zero initial conditions, where F is a smooth function. By taking formally 
the Fourier transform in the x variable, and using the fact that J- C = — \P, 
we obtain that Tu satisfies the following equation: 

d 2 (Tu) 

-fP-tt, = -¥(£).F«(f , + 0, (37) 

with zero initial conditions. Equation ( 137)) can be solved using Duhamel's 
principle. We obtain: 

Fu(t, = -yL= f Bin( - s))TF(s, £)ds. 
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We apply formally the Fourier inversion formula. Multiplying by tp G iS(IR d ), 
and integrating dx, we arrive to the following (formal) definition of a weak 
solution of f )36|) : 

u { t, V )='ff ^§^£% (^F( S , (38) 

If instead of the smooth function F we consider the random noise W, the 
integral above is replaced by a stochastic integral FW(ds, d£), where FW is 
a Gaussian process which we define below. 

As in [2], we let be the completion of with respect to the 

inner product: 



(<Pi,<P2)v(Rd) ■= <Pi(x)<P2(y)f(x -y)dxdy 



^i(£)JV 2 (£Md£). 

Let V(M. d ) be the completion of iS(lR d ) with respect to the inner product: 
(V'l, -02)^5) := {Fip^F^vm- (39) 
Note that if the noise W is white in space, then by Plancherel theorem, 



= (27r) d (^,^ 2 ) i2(Rd) and V(R d ) = V(R d ) = L 



The following lemma gives a more direct way of calculating (ipi, 1P2) 
Lemma 4.1 For any ipi,ip2 G S(M. d ), 

Proof: Note that for any cp G L 1 ^), J>(0 = where 
7-V(0 := / e i$ >(i)di, V£ G K d . 



We denote (p { := F^i G S(R d ) for z = 1, 2. We obtain: 



W>iMzto2s = / / M%)My)f( x -y)dxdy 



F-\F^){i)F-\F^mii{di). 
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By the Fourier inversion theorem, T~ X {T^^ = (27r)% for % = 1,2. The 
result follows. □ 

We endow the space £ with the inner product: 

(hi,h 2 )^p := (J r h 1 ,J r h 2 )np, 

where T denotes the Fourier transform in the x variable. 
Note that by Lemma [4. 1^ for any hi, h 2 G £, 

OO POD 



(h 1 ,h 2 }^ p = a H (2ir) 2d / / / |r - s\ 2H - 2 h x {r, £)h 2 (s, Z)drdsfA(d£). 

JR d JO JO 

For any h G £, we define: 

FW{h) := W{Fh). 
By the isometry property of W, for any h G £, 

E\rw{h)\* = e\w{fk)\ 2 = \\rh\&p = \\h\\y 

Let HV be the completion of £ with respect to the inner product (•, -)^p- 
The map £ 3 h i— )■ JFW{h) G L 2 (fi) is an isometry which can be extended to 
HP. We denote this extension by: 

/■OO 

h i — y 



POO P 

/ / h(t,£)FW(dt,d£)=:FW(h). 

Jo JR d 



This gives the rigorous construction of the isonormal Gaussian process 
J-W = h G ~HV} which was mentioned formally above. 

The following result gives a criterion for a function h to be in HP. 

Lemma 4.2 Let h : M. + x M. d — > C be a deterministic function such that 
h(t, ■) = if t > T . Suppose that h satisfies the following conditions: 

(i) h(t, •) G L 2 (R d ) for all t G [0, T); 

(ii) J-'h G T-LV, where J-'h denotes the Fourier transform in the x variable. 
Then h G H.V and 



a H (2Ti) 2d [ [ [ \r - s\ 2H - 2 h{r,i)h{s^)drds^{di). 

J~R d Jo JO 



nv 



In particular, the stochastic integral of h with respect to the noise TW is 
well-defined. 
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Proof: Let g = J-h. By the Fourier inversion formula on L 2 (R d ), h(s,£) 
(2Tr)~ d J r (s,^), and hence, 



Fg(s,0 = (2«rh(s,Z). (40) 

Since g G HP, there exists a sequence (g n )n>i of the form g n (s, x) = n (s)7„(x), 
where <p n is a linear combination of indicator functions l[o )0 ],o G [0, T] and 
7„ G <S(M d ), such that ||g n - g>||^p -» (see 

Let ^ n := (27r)- d 77^ G S(M d ) and fr„(s,£) = n (»V>n(O- Then 



^n(s, = 0„(s)^ 7 n(O = (27T) d n ( S )^„(O = (27r) d /l n ( S , f). (41) 

Using (J40]) and we obtain that 

a H (2n) 2d [ [ T f \r - s\ 2R - 2 {h n - h)(r,£)(h n - h)(s,Z)drd Sf x(dZ) 
Jm d Jo Jo 

= a H / / |r - s| 2H ~ 2 (Js( n - J r g)(r,£)(J r g n - J r g)(s,£)drdsn(d£) 
Jm. d Jo Jo 

= hn - g\\u V -» 0. 

The conclusion follows. □ 



We now return to equation (|2J). By analogy with (J38J), we say that the 
process {u(t, <p); t > 0, tp G 5(M d )} defined by: 



^ = 1 f'f 3in(^||- S )) ^ 
W 1 Jo V V*(0 

is a u>ea/c solution of (j2J). The stochastic integral above is well-defined if and 
only if the integrand 



1 1 . sin(v^(0(t- S ))- 



(s,0 ^ Kvi s ,£) = Jo~XdMo,t]( s ) TfTTV FcpiZ) belongs to HP. 

To check that G "HP, it suffices to show that h t , v satisfies conditions 

(i) and (ii) of Lemma I4.2I Condition (i) holds since | sinx| < |x| for any x. For 

(ii) , we have to show that g t)ip : = J-h t)l p G HP. For this we apply Theorem 2.1 
of [I]. Note that the function (s, £) >->■ J~9t,<p{s,£) = (27r) d ht >( p(s,£) satisfies 
conditions (i)-(iii) of this theorem. So, if suffices to show that: 

r poo /»oo 

I t>v := a H (2ix) 2d / / \r - s\ 2H - 2 h t ^{r,i)h t ^{s, i)drds^{di) < oo. 

JR d Jo Jo 
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Note that 

h v = I N t (0\^(0\ 2 fi(dO, 

where 

Nt{0 = W) 11 MrVm)M*VWW - s\ 2H - 2 drds. 
Using the fact that | sinx| < \x\ for any x, it follows that 

iVt(0 < a H f frs\r - s\ 2H - 2 drds < t 2H+2 , 
Jo Jo 



and hence I ttV < t 2H+2 J Rd |.F<£>(£)| 2 /^(<i£) < oo. This proves that u(t,(p) 
is a well-defined random variable in L 2 (Q), for any t > and (p £ 
Moreover, 

£Kt,<^)| 2 = £|jw(/^)| 2 = H^ll 2 ^ = / t>v . 

To define the random field solution of ([2]), we proceed as in the case of 
the parabolic equation. We define the norms: 

|M| 2 := E\u(t^)\ 2 = [ N t (0\^(0\ 2 KdO 

r ( 1 \ H+1/2 

Let Mj and Z t be the completions of <S(R ) with respect to the norms || • 
respectively ||| • \\ t . Let M = n t> oM t and Z = n t> oZ t . By Lemma [A. II 
(Appendix A), Z t = Z s = Z for any s, t > 0. 

We say that equation ([2]) has a random field solution if d x E M for any 
x £ M d . In this case, the random field solution is defined by {u(t,x) = 
u(t,5 x );t > 0,x £ M d }. 

The following result gives some upper and lower bounds for N t (^). For 
the upper bound, we use an argument similar to Proposition 3.7 of [I]. For 
the lower bound, we use a new argument. 

Theorem 4.3 For any t > and £ £ R. d , 

H+l/2 / \ H+l/2 



where D^' and are some positive constants depending only on H . 
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Proof: We first prove the upper bound. Suppose that t 2 ^/(^) < 1. Using 
(TO)]) , the fact that IM^i/hyo^ — t 2H _1 IMIi 2 (o t ) and | sinx| < x for all x > 0, 
we obtain: 

1 1 / 1 \^+V2 

- _ft 2 f 2 #+2 < -h 2 2 H+1/2 t I 1 

where for the last inequality, we used the fact that y < f _ 2 ^^ if ^ 2 ^(0 < 1- 
Suppose next that t 2 \I/ (^) > 1. We denote 



Using the change of variable r' = J"v/w(£) and s' = s-\/\I/(£), we obtain: 

^(0 = ^7ml|sm(-)llW)- 

We now use Lemma B.l of [I] for expressing the ?{(0,T)-norm of the 
sinus function in the spectral domain. We obtain that: 

r I |-(2ff-l) 

N *® = *(#m I [/l(T) + ^ (T)]rfT ' (42) 

where 

/t( t ) = sin(rT) — r sinT and gr( r ) = cos(rT) — cosT. 
Letting <p{x) = sinx, we have: (see the proof of Lemma B.l of [4]) 

We split the integral in f|42|) into the regions |r| < 1/2 and |r| > 1/2, and 
denote the two integrals by iV t (£) and iV t (£). Using the same argument 
as in the proof of Proposition 3.7 of [4|, we get: 

iV^(0 < C , C f„ - ■ < C^-1 2U -H ( ) , (43) 

* ~ ^(o H+1 i-H~ i-H \i/t 2 + *(0y 
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where C = 11.11, and for the second inequality we used the fact that ^7^172 < 
^and^ ift 2 *(0>l- 



On the other hand, 



where 



fUr) + 9 2 t(t) 



I(T):= / J1 ) ' ' dr = 2n / sin 2 xdx = ttT 

(t 2 — lr 



1 



sin(2T) 



(44) 



2T 

by Plancherel's theorem. This yields the estimate I(T) < 2nT. We obtain 



Ni 2) {t)<c H tf H - 1 2%t(-±-r) <c n 2™- x 2id 



H+l/2 / o \ H+l/2 



*(0/ " vV* 2 + *(0 

where for the second inequality we used the fact that t 2, 3/(£) > 1. 
Combining (T43]) and ( 1451) . we conclude that: 



(45) 



H+l/2 



The upper bound follows, letting 



= max 



fl fe 2 2 ff + l/2 cr_£fl_ 2 3H-l/2\ 

\3 H ' 1-H J 



We now treat the lower bound. Suppose first that t 2 \I / (^) < 1. Using the 
fact that sinx > xsin 1 for all x G [0, 1], we obtain: 



N t (0 > a H sm 2 l [ [ rs\r - s\ 2H ~ 2 drds = a H sm 2 l ^lE ^ +2 
Jo Jo H + 1 

. 2 J(2,2H-1) J 1 V 
> a«- sin 1 1 — — 



H+l/2 



where j5 denotes the Beta function and we used the fact that t 2 > pr^^y. 

Suppose next that T 2 = t 2 ^/(^) > 1. Let p > 1 be a constant which will 
be specified below. We use f|42|) . Since the integrand is non-negative, N t (£) 
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is bounded below by the integral over the region |r| < p. In that region, 



m - l /(T) " L_ P (r*-!)' dr ) • (46) 

Using (jHJ) and the inequality 1 — (sina;)/a; > 1/2 for any a; > 2, we get 
the lower bound: 

I{T)>-T, since T > 1. (47) 

To find an upper bound for the second integral in the right-hand side of 
f )46|) . we use the fact that: 

/|(r)+^(r)<2T(l + |r|) 2 , VrGl. 

(To see this, note that |/t(t)| < 1 + \ T \ an d |/t(t)| < 2T|r|, since 
| sin a; | < \x\. Hence, fxij) < 2T|r|(l + \t\). Similarly, |g:r(r)| < 2 and 
|#r( r )l < T(l + |r|), since |1 - cosx| < |ar|. Hence, g^(r) < 2T(1 + |r|).) 



It follows that: 

r|>p (^-l) 2 



dr < C P T, (48) 



where C p = 2 Jj r|>/9 dr. Using flU, (STJ) and gHD, we obtain that: 

^«)^%^(I-^)'V*K). 

Choose p large enough such that C p < ir/2, e.g. p = 4, for which C p < 4/3. 
Using the fact that > t -2^^ , we get 

iVt(e) > c g 4-( 2g - 1 ) f- - -V ( - - 1 

- V 2 3y V 1 /* 



'(0 



The lower bound follows, letting 



□ 

A consequence of the previous result is that the norms || ■ \\ t and 
are equivalent, for any t > 0. 



k 2 
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Corollary 4.4 For any t>0,<pe S(R d ), 

d H t£{t 2 ;(p) < E\u{t,if)\ 2 < D H t£(t 2 ;<f). 

Hence M t = Z t i for any t > 0, and M = Z. 

Below is the main result of this section. 

Theorem 4.5 (Existence of Solution in the Hyperbolic Case) Assume that 
(G?3J) holds. Then has a random field solution if and only if holds. 

Proof: As in the proof of Theorem 13. 3[ one can show that ([H]) is a necessary 
and sufficient condition for 5 X G Z\ = M for all x e M. d . We omit the details. 
□. 

Example 4.6 (Stable processes) As in Example 13.51 let C = —(—A)' 3 / 2 for 
P e (0,2]. Then = c^f and ® becomes: 

r ( i \ H+1/2 

When f(x) = c a) ri|a;|~( d ~ a ) with < a < d, ( )49l) is equivalent to 

(^H + ^ J p > d - a, 

whereas for f{x) = Ylf=i( a Hi \xi\ 2Hi ~ 2 ), (SHJ) is equivalent to 



Remark 4.7 (Fractional Powers of the Laplacian) As in Remark I3.6[ we 
can consider the case C = —(—A)' 3 / 2 for arbitrary (3 > 0. Note that the 
fundamental solution G of duu — Cu = exists and satisfies: 



(see p. 11 of [H]). Using Theorem 2.1 of [I] and estimates similar to those 
given by Theorem 14.31 above, one can show that a random field solution of 
02]) (in the sense of [I]) exists if and only if f T4"9]) holds. 
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A Some elementary inequalities 

Lemma A.l For any a, (3 > 0, 

Cl (a,P) 2H T* H (P) < T* H (a) < c 2 (a, (3) 2H T* H (f3), 

where ci(a,/3) = (/3 A l)/(a + 1) and c 2 (a, (3) = (/? + l)[(l/a) V 1]. 

Proof: We denote by T* Hl (a) and T* H2 {ot) the integrals over the regions 
{2Re#(£) < 1}, respectively {2Retf(f)' > !}■ Using the inequality (a + 
I)- 1 < [1 + 2Re#(0r 1 < a" 1 if 2Retf (f) < 1, we obtain that: 

Combining this with the similar inequality for T* H1 (f3), we get: 

(^i) 2 t h,M < ^hM < {^Y^kM- (50) 

A similar argument works for T* H2 {ce). We obtain: 
/ 1 \ 2H 

T ^ 2(/3) - T * H ' 2(a) ~ {(5 + 1) T ^(/3)- (51) 

The result follows by taking the sum of ( l50l) and (j5Tjl . □ 

We recall the definitions of the functionals S(t; tp) introduced in Section 
[31 respectively Section HI To make a distinction between these functionals in 
the two cases, we use the index "par" for parabolic, and "hyp" for hyperbolic: 

£ ^ = I( iA + Rew(o ) g+ ' /2 '^K)lV(^). 

Lemma A. 2 For an?/ s > 0,t > and <£> G iS(IR d ) ; 

Ci(s,t) 2H £ par (,s;y?) < £ par (t;y?) < c 2 (s,t) 2H £ par (,s; ip) 
Cl (s,t) H+1 / 2 S hyp (s;p) < S hyp (t;<p) < c 2 (s,t) H+1 / 2 £ hyp (s;ip), 
where c x (s,t) = (s^ 1 A l)/^" 1 + 1) and c 2 {t) = (s' 1 + l)(t V 1). 
Proof: The argument is similar to the proof of Lemma IA.1I □ 
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B The 7^(0, T)-norm of the exponential 

The next result gives the expression of the "H(0,T)-norm of the complex- 
valued exponential function in the spectral domain, which is needed in the 
proof of Theorem 13.11 

Lemma B.l Let (p(x) = e - x ( a+lb ) for x G (0,T), where a, b G R. Then 



where c H = T(2H + 1) sm(nH)/(2H). 

Proof: For complex- valued functions (p G L 2 (0, t), we can apply Lemma A.l 
of jl] to Keep and Imip to obtain that: 



where frij) = sin[(r + b)T] and grij) = e at — cos[(r + b)T]. The result 
follows by (|52jl. □ 

C A version of Plancherel theorem 

The following result is a version of Plancherel theorem needed for the calcu- 
lation of (P r+S (f * <f>))(x) in the proof of Lemma [3.101 

Lemma C.l For any <p G L 2 (R d ), fa G L 2 (R d ) and fa G L 1 (R d ) n L 2 (R d ), 

I I fa(x)fa(y)<p(x - y)dydx = — ^ / FfaffiTifolg) F<p{£)d£. 

Proof: By Young's inequality, \\<p * fa\\ 2 < IMhll^Hi and (f * t/> 2 G L 2 (IR d ). 
The result follows by Plancherel theorem, since 





(52) 



where Fq,t<p{t) := J e lTX ip(x)dx. 

An elementary calculation shows that for ip{x) = e ~ x ( a + lb )^ 




□ 
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